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ABSTRACT 

Cyclotron resonance scattering features observed in the spectra of some X-ray pulsars show a significant 
changes of the line energy with the pulsar luminosity. At high luminosities, these variations are often associated 
with the onset and growth of the accretion column, which is believed to be the origin of the observed emission 
and of the cyclotron lines. However, this scenario inevitably implies large gradient of the magnetic field strength 
within the line-forming region, which makes the formation of the observed line-like features problematic. 
Moreover, the observed variation of the cyclotron line energy is much smaller than could be anticipated for the 
corresponding luminosity changes. We argue here that a more physically realistic situation is that the cyclotron 
line forms when the radiation emitted by the accretion column is reflected from the neutron star surface, where 
the gradient of the magnetic field strength is significantly smaller. We develop here the reflection model and 
apply it to explain the observed variations of the cyclotron line energy in a bright X-ray pulsar V 0332+53 over 
a wide range of luminosities. 

pulsars: general - relativistic processes - scattering - stars: neutron - 



Subject headings: line: formation - 
X-rays: binaries 

1. INTRODUCTION 

X-ray pulsars are neutron stars in binary systems accret- 
ing matter usually from a massive companion. These neutron 
stars have a sufficiently strong magnetic field, which channels 
accreting gas towards magnetic poles. Strong magnetic field 
modifies the observed X-ray spectrum often manifesting as 
the line-like absorption features, the so-called cyclotron lines. 
Such cyclotron resonance scattering features (CRSF), some- 
times also with harm onics, are observed in the spectra of sev- 
eral X-ray pulsars (ICoburn et al.l l2002; Filip pova et al.ll2005t 
ICaballero & Wilmsll2012T) ~ 

In some cases, the luminosity related changes of the line 
energy are observed, suggesting that configuration of the line- 
forming region depends on the accretion rate. The line en- 
ergy has been reported to be positively (in relatively low- 
luminosity sources; see Staubert et al. 2007; Klochkov et al. 
120121) and negatively-c orrelated with luminosity (in high- 
lumin osity sources; see Miha ra et al.| [l998; Tsygank ov et al.l 
2006), as well as uncorrected with it ((Caballero et al. 2013). 

This diversity is yet to be explained, and in this work we 
will focus only on the high-luminosity case. The negative 
correlation of the CRSF energy with luminosity here is usu- 
ally explained with the ons et and growth of the ac cretion col- 
umn at high luminosities (Basko & Sunvaevlll976l) . In this 
scenario, the height of the column, and, therefore, the aver- 
age displacement of the emission and the line-forming regions 
from the neutron star surface increase with luminosity. The 
magnetic field weakens rapidly with distance from the neu- 
tron star and, therefore, the CRSF should shift to lower ener- 
gies. The problem is, however, that the predicted shift is much 
larger than the observed on e. The column height dep ends on 
luminosity almost linearly (Basko & Sunyaev 1976) and the 
magnetic field weakens with distance as r~ 3 , and yet bright- 



ening by more than an order of magnitude yields at most 25% 
decrease in the CRSF energy (Tsygankov et al. 2006, 2010). 
Moreover, large gradient of the magnetic field and of the ac- 
cretion velocity are expected to smear out the line-like fea- 
tures making it difficult to explain why we observe CRSFs at 
all. 

Sev eral authors (Burderi et al. 2000; Krevkenboh m et ail 
2008) considered variation of the magnetic field along the 
neutron star surface as a possibility to explain the observed 
variation of the cyclotron line energy with the pulse phase, 
which has similar magnitude as those associated with the lu- 
minosity changes. Interaction of accreting plasma with the 
magnetosphere of the neutron star defines the geometry of 
the emission region, so one could imagine that a change in 
the accretion rate could also offset the location of the polar 
cap, leading to a change in the observed CRSF energy. How- 
ever, the observed variations would imply an unrealistically 
large shift of the hotspot location (by 50-60 degrees) from 
the magnetic pole. Furthermore, the observed change in the 
CRSF energy in this case would imply a dramatic change in 
the pulse profile shape, which is not observed. 

On the other hand, a significant part of the column ra- 
diation should be intercepted by t he stellar surface be- 
cause of the r elativistic beaming (iKaminker et al.l 119761: 
Lyubarskii & Sunyaev 1988). Small variations of the i3-field 
strength over the surface (only by a factor of two in the case of 
the dipole field) would imply that the variations of the typical 
cyclotron energy are small. Thus, it is more natural to assume 
that the line is formed in the atmosphere of the neutron star 
due to reflection of the intercepted radiation. Increase in the 
mass accretion rate results in a larger luminosity, a higher col- 
umn height and a larger illuminated part of the stellar surface 
(see Fig.[TJ. If the B-field decreases away from the magnetic 
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Figure 1. Sketch of the accreting X-ray pulsar geometry and the emergent spectrum. The larger is the accretion rate, the higher is the column, the larger 
illuminated fraction of the stellar surface, the weaker the average magnetic field, and the smaller the cyclotron line energy. 



poles, the cyclotron line energy should then decrease with the 
luminosity and a negative correlation between the luminosity 
and the cyclotron line energy is reproduced. Here we dis- 
cuss this scenario quantitatively and compare the model pre- 
dictions with the data taking during a bright outburst of the 
transient X-ray pulsar V 0332+53. 

2. MODEL SET UP 

Let us start from the physical picture of the accretion 
on the magnetized neutron star follow ing ear lier papers by 
iBasko & Sunyaevl (|1976|) and Kamink eretai] (QI976). It is 
possible to distinguish the two regimes of accretion onto 
magnetized neutron stars depending on the mass accretion 
rates. At low accretion rate, free-falling protons heat part 
of the neutron star surface near its magnetic poles, and these 
bright spots radiate energy in the X-ray range. At high ac- 
cretion rate, radiation pressure becomes significant and stops 
the infalling material above the neutron star surface in the 
radiation-dominated shock. Below the shock, the matter 
slowly sinks down as the excess emission supporting the col- 
umn escapes through the side walls. The column is expected 
to arise as soon as the luminosity exceeds a critical value 
(IBasko & Sunyaevlll976l) : 
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where n\\ is the electron scattering opacity along the magnetic 
field, kt is the Thomson opacity, M and R are the mass and 
the radius of the star, I is the length of the accretion arc at 
the stellar surface. It is worth noticing that the optical depth 
across the column is of the order L/L* and L* is much larger 
than the Eddington luminosity scaled to the area of the foot- 
print of the accretion column and depends on the accretion 
flow geometry. 

The column height depends on the accretio n rate 
(IBasko & Sunyaevl 1976t lLyubarskii & Sunyaevl 1988b : 
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d is the thickness of the accretion arc, and m = L/L** is a 
ratio of the X-ray pulsar luminosity to the limiting luminosity 
for the magnetized neutron star 
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which corresponds to the column height of h > R. The 
height, where matter stops, varies inside the accretion channel 
and depends on the distance from its borders because the radi- 
ation energy density drops off sharply towards the edge of the 
column. The height has its maximum value near the middle of 
the channel and decreases towards the borders. Therefore, the 
radiation from the already stopped matter should pass through 
a layer of the rapidly falling plasma. These outer layers of 
the column are not supported by the radiation and fall with 
velocity close to the free-fall velocity (3 = v/c — yVs/V 
(here r§ = 2GM/c 2 is the Schwarzschild radius). The op- 
tical thickness of these layers is high enough to change sig- 
nificantly the angular distribution of the emergent radiation. 
As a result, the radiation is directed mainly towards the stel- 
lar surface due to the relativistic beaming. For the electron- 
scattering dominated column, the angular distribution of the 
column luminosity in the laboratory frame i s given by (see 
Appe ndixlAl and Kaminkeret al. 1976; Mitro fanov & Tsyganl 
1978): 
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where a is the angle betwee n the ph oton momentum and 
the velocity vector, 7 = 1/yl — f3 2 is the Lorentz factor, 
D = 1/[7(1 — (3 cos a)] is the Doppler factor, and I is the 
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Figure 2. Fraction of the captured radiation from a point source above the 
neutron star surface as a function of the height-to-radius ratio. The dashed 
blue curve corresponds to the isotropic source in fiat space-time, the dot- 
ted red curve is for the isotropic source accounting for light bending in 
Schwarzschild metric, and the solid black curve is for the emission pattern 
given by equation {5} and accounting for light bending. Here R = 3rs ■ 
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Figure 3. Dependence of the normalized flux on co-latitude for an accretion 
column emitting according to the law l(5] with constant luminosity per unit 
height (i.e. constant Iq). The flux is normalized as f F (9) sin 9 d9 = 1. 
The red solid lines correspond to Ah/h = 1 and the blue dashed lines to 
Ah/h = 0. Different lines show the distribution for various column heights: 
h/R = 0.5, 1, 2 and 4 (from left to right). Here R = 3r s . 



normalization constant. The emission pattern becomes more 
isotropic at large r, because of the radial dependence of the 
free-fall velocity. 

The energy release and the distribution of the emerging ra- 
diation flux along the accretion column depends on the details 
of radiation-matter interaction and is a complicated radiation- 
hydrodynamical problem. For simplicity, we approximate the 
geometry of the accretion column by a thin stick on the mag- 
netic pole. This approximation is reasonably accurate for high 
B-field pulsars up to h < R. Here we also assume that most 
energy is emitted in a region of characteristic scale Ah situ- 
ated above the neutron star surface at height h (see Fig. [TJ, 
so Ah/h = 1 for a homogeneously emitting column, or 
Ah/h = if all the energy is emitted within a thin shock 
at the top of the column. 

Results weakly depend on the compactness of the star, 
which we fix here at R/rs — 3 (corresponding to the radius 
of 13 km for a 1.5M Q neutron star, consistent wit h the recent 
measurements from X-ray bursters, see e.g. Sulei manov et al.1 
l20llh . 

3. REFLECTED FRACTION AND SURFACE FLUX 
DISTRIBUTION 

At a high accretion rate, the accretion column radiates the 
released gravitational energy of the falling matter through its 
sides. Part of the radiation leaves the system directly and part 
of it is captured by the stellar surface. The fraction of the cap- 
tured radiation L c /L depends on the height of the column, 
brightness distribution over the column, radiation beam pat- 
tern and the compactness of the star. In Schwarzschild metric, 
for an isotropic, point-like source at radius r, it is given by a 
simple formula (see AppendixlBl: 
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In case of flat space-time, we can substitute rs = 0. The 
captured fraction is slightly higher when light bending is ac- 
counted for. It drops from the maximal value of 0.5 to 0.1 
when the source rises to h ~ R (see Fig . [21) . 

For the emission pattern given by equation (0 and used 
in the following calculations, the captured fraction increases 
even further, because of the strong beaming towards the star 
(see top curve in Fig. |2). For a higher source position, L c /L 
drops below the isotropic case at h ~ 2R, because of a rather 
small beaming resulting from a smaller free-fall velocity and 
a factor sin a, which reduces the radiation directed along the 
column. 

Variations in the emission height lead also to a dramatic 
change in the distribution of the captured flux F{9) over the 
stellar surface. We illustrate this effect here by considering 
a homogeneously emitting column (i.e. Iq = const in equa- 
tion[5]) of various heights and Ah/h (see Fig. [3]). The details 
of the calculations can be found in AppendixlBl For a low col- 
umn, most of the emission hits the surface in the direct vicin- 
ity of the column at a co-latitude 9 < y/2h/R and the flux 
drops rapidly with distance from the column. For Ah/h = 1, 
the flux F(9) diverges as 1/9, which is a consequence of our 
assumption of an infinitely thin column. When h sw R, the 
star is irradiated almost up to the equator and for even higher 
column, most of the stellar surface receives some flux. In a 
realistic pulsar, the captured flux is a sum of the contribution 
from the two antipodal columns. 

4. FORMATION OF THE CYCLOTRON LINE 

Radiation from the accretion column hits the neutron star 
surface and is reflected. The most important process affecting 
the spectrum of the reflected radiation is Compton scattering. 
Under the conditions of the neutron star atmosphere, the den- 
sity is too small to produce any true absorption of the hard 
X-ray photons either in the continuum or at the cyclotron line 
energy. The cross-section for Compton scattering in strong 
magnetic field is energy-dependent and has strong resonances 
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where b = B/B CI is the i?-field strength in units of the critical 
field strength B„ = mlc 3 /eh = 4.412 x 10 13 G, £ is the 
angle between the field and photon momentum, and m c is 
the electron mass. For ( < 1, the resonance energies are 
E/cs /m c c 2 s=s bn. 

Photons at the resonance energies cannot penetrate deep 
into the atmosphere, they interact in the surface layers and 
scatter back. In the cyclotron line wings, photons penetrate 
deeper into the atmosphere and scatter there with some en- 
ergy shift. If they scatter into the resonance energy, they 
cannot leave the atmosphere because of the larger optical 
depth there and escape instead in the line wings. Thus, 
the lack of the photons near the resonance is not filled in. 
The absorption feature at the resonance energy and the emis- 
sion features in the wings appear in the spectrum of the re- 
flected radiation. The energy separation between the emis- 
sion peaks can rea ch many Doppler widths (see Avrett 1965, 
lAvrett & Hummenll965l and Sect. 6.4 of lIvanovlll973i) . Be- 
cause the emission lines are broad, they merge with the con- 
tinuum and cannot be easily separated. For resonant Comp- 
ton scattering in magnetic field, the line is not symmetric and 
the red wing of the emission line is stronger, because of the 
recoil ( Wasserman & Salpeter 1980; Arava & Harding 1999; 
Araya-Gochez & Harding 2000; Harding & Lai 2006), which 
gives the relative photon energy shift of about AE/E ~ 
—E les /m a c 2 ps —b. The absorption features at the harmonic 
energies can be even stronger, because the absorbed photons 
are mostly reemitted at the energy of the fundamental. The 
typical X-ray pulsar spectrum cuts off at ^30-50 keV and 
therefore the contribution of the harmonics will be negligible 
in high-field pulsars (as is the case of V 0332+53). 

Neglecting the asymmetry in the line shape, the centroid of 
the CRSF in the reflected spectrum averaged over the surface 
and all angles is determined by the B-field strength weighted 
with the distribution of flux, F(9) (Pig.[3j, and the line equiv- 
alent width, EW(0), over the neutron star surface: 
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where z is the gravitational redshift. Assuming that the line 
EW in the reflected radiation is constant over the surface, we 
can relate the cyclotron line centroid to the mean field as 
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Variations of the magnetic field over the surface lead to the 
smearing of the line and its minimum width is then related to 



If (a b /(B)) x E cyc \ is smaller or comparable to the separation 
between the emission peaks, the CRSF will remain strong in 
the total spectrum of the reflected radiation. 

At low accretion rate, the column is very low, h/R -C 1, 
and only the polar region is illuminated, resulting in the line 
energy E cyc \ — Eq corresponding to the polar field Bq. Be- 
cause for the dipole field, the S-field strength drops only by a 
factor of two from the pole to the equator, there is an obvious 
limitation for changes in the cyclotron line energy: the line 
should lie in the energy interval [Eq/2; Eq]. A more realis- 
tic lower limit on the line energy can be obtained assuming 
a uniformly illuminated surface, i.e. F(ff) = const. In that 



7r 

(B) min = jB J sin 0\/l + 3cos 2 0d0 ~ 0.7B . (13) 

o 

Therefore, the line should not change its energy significantly 
and the model predicts the range for the l ine centroids of 
\Q.7En\Ep\ in agreement with observations (Tsygank ovet al.l 
l2Q06ll20Toh . 

In Fig. IH we present theoretical dependences of the mean 
magnetic field, standard deviation, and the captured fraction 
on the height of the column for different values of Ah/h. 
These quantities will have a similar behavior as a function 
of luminosity too. We see that for the point-like source, vari- 
ations in all quantities are the largest, while for a homoge- 
neously emitting column, they are the smallest. If the whole 
column emits (i.e. Ah/h = 1), the polar regions are al- 
ways strongly illuminated and (B) varies by only 10%. On 
the other hand, for Ah/h = 0, the mean field varies by at 
least 25%. A negative correlation between the column height 
(luminosity) and the cyclotron line energy is produced. The 
smearing of the line, because of the magnetic field varia- 
tion, does not reach more that 10-15% (see Fig. I4J5). Thus 
if the emission peaks in the reflected radiation are separated 
by more than 20%, the CRSF would remain strong even at 
high h/R. The strength of the CRSF in the total spectrum de- 
pends also linearly on the captured luminosityfraction and is 
expected to decrease at high column (see Fig.0};); therefore, a 
correlation between the line equivalent width and the energy 
is expected. 

5. COMPARISON WITH OBSERVATIONS 

We would like to compare our model with the data 
from the X-ray pulsar V 0332+53 obtained with the RXTE 
and INTEGRAL observatories during outburst in 2004-2005 
(Tsyg ankov et al.l [20061 120 10). This source shows negative 
correlation between the luminosity and the energy of the cy- 
clotron line. It is the only data set which has the information 
about the behavior of the object in a so wide range of the lu- 
minosities, from ~ 10 37 up to ^ 4 x 10 38 erg s _1 . 

For a given parameter of the accretion column Ah/h, we 
compute the distribution of the reflected flux F(9) as a func- 
tion of the column height h. This is then converted to the 
average magnetic field (B) in units of the polar one Bq us- 
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Figure 4. Theoretical dependences of (a) the mean magnetic field in the 
refle cting zone for dipole field (in units of polar field) as defined by equation 
tilt , (b) the magnetic field standard deviation given by equation 1121 . and (c) 
the fraction of the captured luminosity on the h/R, ratio for R = 3rg. Black 
solid, red dotted, blue dashed, green dot-dashed and pink triple-dot-dashed 
line correspond to Ah/h=0, 0.3, 0.5, 0.7, 1.0, respectively. 
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Figure 5. Dependence of the cyclotron line energy on the l uminosity in the 
X-ray pulsar V 0332+53 (circles with the error bars; from Tsv gan kov et all 
2010) and the best-fit theoretical relation (solid line) for Ah/h = 0.1, r] = 
15, E = 29.5 keV, and L** = 21 X lO^ 8 erg s _1 . 



Table 1 

Best-fit parameters 



Ah/h 


E 
(keV) 


L** 
(10 38 ergs -1 ) 


x 2a 


0.0 


29.5±0.2 


22.3+f* 


96.4 


0.1 


29.5±0.2 


20.9±1.8 


96.4 


0.2 


29.5±0.2 


19.41^ 


97.3 


0.3 


29.6±0.2 


17.8±1.6 


98.5 


0.4 


29.6±0.2 


16.1±1.6 


100.0 


0.5 


29.6±0.2 


14.3±1.6 


103.7 


0.6 


29.6±0.3 


12.4+1-1 


108.4 


0.7 
0.8 
0.9 
1.0 


79 7+0.5 
zy -'-0.3 
29 9+10 
zyy -0.3 


10.2+J;g 

, ,+0.35 
'•'-0.13 

,2 5+2.4 


116.7 
130.6 
124.2 

227.3 



d For 94 degrees of freedom. 

ing equation ( fTTT i assuming a dipole field and to the cyclotron 
line energy E cyc \ in units of the polar value Eq using equation 
( [Tol l. Transition from the dependence on h to the dependence 
on the luminosity is made as discussed in Sect, fusing equa- 
tion (0, where 77 and L** are the parameters. 

The cyclotron-line energy depe ndence on the lu minosity 
for the X-ray pulsar V 0332+53 dTsv gankov et all 120101) is 
shown in Fig. [5] It is clear that the data at luminosities above 
3.6 x 10 38 erg s _1 cannot be described by our model and we 
neglect them in the fits. The remaining data are fitted with 
two free parameters Eq and L**. We fix ry = 15 (see equation 
[3]), because the results depend on that parameter very weakly 
and vary Ah/h in the range between and 1 with the step 0. 1 . 
Taking the errors on measured E cyc \ at their face values, gives 
the reduced \ 2 significantly above unity, because of a large 
spread of the line energies at a fixed luminosity. Therefore, we 
add a systematic error of 0. 15 keV in quadrature to the statis- 
tical errors of E cyc \. A good agreement with the data, with the 
X 2 values below 100 (for 94 degrees of freedom), is achieved 
for E « 29.5 keV, Ah/h < 0.4 and L"k2x 10 39 erg s" 1 
(see Table [TJ. The \ 2 grows with Ah/h and L** decreases 
from about 2.2 x 10 39 erg s _1 to 10 39 erg s _1 when Ah/h 



Poutanen et al. 



changes from to 0.7. The values of L** obtained in the fits 
should be taken as the upper limits on the actual value given 
by equation (0), because the radiation from the upper part of 
the column might be completely blocked by the falling mate- 
rial and photons escape only at a height significantly smaller 
than the top of the accretion shock. 

The fitting procedure seems to indicate that the best de- 
scription of the data is achieved for a column with most of 
the emission coming from its top. This does not necessarily 
mean that the lower part is not emitting, but rather that this 
emission does not hit the neutron star surface. It is actually 
expected, because in the lower part of the column the optical 
depth of the falling plasma above the shock becomes smaller 
and the photons from the shock freely escape to the observer. 
Furthermore, our results are based on the assumption of the 
constant with latitude equivalent width. In reality, the scatter- 
ing opacity is a strong function of the angle betw een the pho- 
ton momentum and the magn etic field (see e.g. Pavl ov etal] 
Il980t ISuleima nov et al.ll2012l) . which necessarily will result 
in some latitude dependence and possibly will increase the 
line EW of the spectrum reflected from the equatorial region. 
In that case, our result of the top-dominated emission from 
the column might be an artifact of this effect. 

The best-fit parameter indicate that the height of the column 
reaches R at L ~ 4 x 10 38 erg s _1 . This is close to the lu- 
minosity, where large deviations of the model from the data 
are visible. The possible reason for this mismatch is that our 
approximation of the accretion column by a thin stick breaks 
down. In addition, at such large L, radiation from the high 
anti-podal column starts to hit the equatorial region and be- 
comes observable. This reduces the average "observed" B- 
field. In our model, the presented quantities are averaged over 
all observer angles and this effect is not accounted for. 



6. SUMMARY 

In this paper, we have proposed a reflection model for the 
cyclotron line formation in X-ray pulsars. At high accretion 
rates, the accretion column is predicted to have a significant 
height illuminating a large fraction of the neutron star surface. 
Strong beaming of the column radiation towards the surface 
guaranties that a significant fraction of the column emission 
is reflected from the surface. 

We argued that the reflected spectrum should have a strong 
CRSF at the energy depending on the local magnetic field 
strength. Small variations of the magnetic field along the sur- 
face imply that the line centroid energy can vary by at most 
30%, from the value corresponding to the field at the pole to 
the whole surface-averaged field. Changes in the pulsar lu- 
minosity are expected to correlate with the illuminated frac- 
tion of the stellar surface, and anti-correlate with the average 
magnetic field and, therefore, with the line centroid energy, 
exactly as observed during the outburst of the X-ray pulsar 
V 0332+53. Our model has profound implications for the 
interpretation of the data on the cyclotron lines observed in 
X-ray pulsars. 

In order to predict the line parameters more accurately, a 
detailed model of the reflection of the column radiation from 
the atmosphere is required. It necessarily should include the 
dependence on the orientation of the local B-field as well as 
the angles of the incoming and reflected radiation. We would 
then be in a position to predict variations of the line energy 
and EW with the pulsar phase. This project is outside of the 
scope of the present paper and is left for a future study. 
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APPENDIX 

A. EMISSION FROM THE RAPIDLY FALLING PLASMA IN AN ACCRETION COLUMN 

Let us consider an accretion column with the gas falling with velocity f3 = v/c. Let us assume that the angular distribution of 
radiation escaping from the surface in the gas comoving frame is given by the bolometric intensity 



7'(C') = Io(l + acosC / ), 



(Al) 



where £' is the angle measured from the local normal to the surface, Iq is a constant, and a is the anisotropy parameter. In the 
calculations, we have chosen a = 2 , which corresponds to the case of electron-scattering dominated optically thick atmosphere 
(Chandrasekhar 1960; Sobolev 1963), but keep here the formulae as general as possible. We define the coordinate system with 
the z-axis along the direction of motion of the gas, and the x-axis along the normal to the surface. The four- vector of the photon 
momentum in the lab frame is k = k {1, u>}, where k — hv/m c c 2 and u> is the unit vector along the photon momentum, which 
makes angle a with the z-axis: u> = (sin a cos (?!>, sinasin</>, cosa). Relative to the surface normal n = (1, 0,0), it makes angle 
£, so that cos £ = u> ■ n = sin a cos <f>. 

The photon four- momentum in the frame comoving with the spot, k!_ = k' {1, a/}, is obtained from the Lorentz transformation. 
The energy is transformed as k' = k/D and the unit vector along photon momentum is 



u> 



Here the Doppler factor is 



D sin a cos 6, D sin a sin < 



D = 



cos a 



P 



1 — (3 cos a 



(A2) 



(A3) 



7(1 — /3cosa) ' 

and 7 = l/\/l — /3 2 is the Lorentz factor. Note that because of the relativistic aberration, the cosine of the projection angle we 
see the surface element that moves along its surface is (see also Poutanen & Gierlihski 2003; Poutanen & Beloborodov 2006): 

cosC' = Dcos( — D sin a cos (f>. (A4) 
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AT 

^-= /'(C) cos C. (A5) 



The bolometric power emitted by a unit area surface element in the comoving frame per unit solid angle is given by 

dL 

du> 
The angle-integrated emitted power is 

L' = J> ?—du>' = f dip' J J'(C') cosC'dcosC' = /ott (l + |a) ■ (A6) 

In order to evaluate the flux on the neutron star surface, we need to know the emitted power in the lab frame, which is given by 

dRvbicki & Lightmanl 19791) 

dL D 3 dL' D\,,^. . 

— = — = — /'(C) smacos^. (A7) 

au 7 ou' 7 

Note that here the first transformation does n ot involve factor D A , because we consider a steady-state source, not a moving blob 
dRvbicki & Lightmanll97 9: Sikora eTalJ 19971) . For a thin accretion column, photons emitted at any azimuthal angles <fi will have 
basically the same trajectory that is described just by angle a. Thus, we can integrate over to obtain the final expression for the 
emission pattern 

dL r /2 dL D 4 ( ir \ 

- = / dq> —— = Iq — 2 sin a 1 + a — D sin a) . (A8) 

a cos a J—n/2 du j V 4 / 

This expression is different somewhat from that given by Kaminker et al. (1976) and Mitrofanov & Tsygan (1 1978b as it contains 
one less Doppler factor. For the emission forward-back symmetric in the comoving f rame , the total emitted power is Lorentz 
invariant (Rybicki & Lightman 1979), which can be checked by integrating expression ( 



/ A T 
- dcosa = L'. (A9) 

_! a cos a 

B. PHOTON PROPAGATION IN THE VICINITY OF THE NEUTRON STAR 

Figure|6]depicts the trajectory of a photon (red solid line) which is emitted at the accretion column at height h above the neutron 
star surface in the direction that makes angle a with the radial direction towards the stellar center. The photon's path, described 
by the distance r = R + h and azimuthal angle ip, obeys the equation of motion (e.g. Misne r et al.lll973b iPechenick et al. 1983): 

%) +( i -«),,> = i. pi) 

where u — rs/r and b is the impact parameter in units of the Schwarzschild radius rs = 2GM/c 2 for a neutron st ar with mass M. 
The impact parameter and the angle, a, between the radial direction and the photon trajectory are related by (e.g. Pec henick et al.l 
1983; Beloborodov 2002) ' sinQ , 

b= 7= (B2) 

uyl —u 
The azimuthal angl e ip measured from the continuation of the trajectory towards infinity (dashed red line in Fig.[6]l can be obtained 

from equation (IB lb : 

V>(r, a) = J du' [b- 2 - (1 - u')u' 2 } ~ 1/2 . (B3) 

For r larger than 2rs the e lliptic integral ( IB3t can be approximated with a high accuracy by (Belobor odovl 120021: 
iPoutanen & Beloborodovl2006l) 

cos a — u(r) + [1 — u(r)] cos ip. (B4) 

The photon trajector y wi th impact parameter b meets the neutron star surface of radius R at an angle a.Q relative to the normal, 
which is (see equation (IB2I )): 

u{R)Jl-u{R) 

sinao = sin a — (B5) 

u(r)yl — u(r) 

The intersection point lies at the azim uth t />q that can be obtained from equation dB4b substituting there compactness u = u(R) = 
rg/R and «o instead of a. Equation dB5l ) immediately constrains the maximum angle a, when the trajectory is still intercepted 
by the star, corresponding to a$ = tt/2: 

R I 1 - u(r) 

sma max = — a — -. (B6) 

r y 1 — u(R) 

For a source at distance r from the star emitting according to the angular pattern dL(a)/d cos a with total luminosity L, we then 
can obtain the fraction of radiation captured by the star: 

-^C J COS (Umax U COS Ci /"D'7\ 

-£- = — ^ — jrj^\ . • wo 



I 1 

J COS Oiu 


dL(a) , 

, - ; a cos a 

lax d cos a 


/i 


dL(a) , 
dcosa dc0Sa 
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Figure 6. Geometry of light bending in Schwarzschild metric. The red solid line shows the photon trajectory from the emission point at distance r 
the stellar surface at radius R and co-latitude 9. The red dashed line is the continuation of the trajectory in the opposite direction to the infinity. 



For an isotropic source, equation (IB7) is reduced to 

Le = 1 

L ~ 2 
Distribution of the intercepted luminosity over the neutron star surface is given by 



(1 



cos a r 



0- 



R + hto 



(B8) 



dL c (6) dL(a) dcosa dL c 

~dL 



■ n , in if * (B9) 

dcosu dcosa dcosu dL 

where the first factor on the rhs is the specified emission pattern. We compute the second factor num erically differentiating 

9 = tpo — if) as a function of a, with both if> and ipo obtained via approximate light bending equation (IB4b and a^ fro m equ ation 

( IB5b - The last factor, dL c /dL = [1 — u(r)]/[l — u(R)], just contains two redshift factors of the type 1 + z = 1/Vl — u: one 

from the photon energy change in the gravitational field and another one from the time dilation. 



REFERENCES 



Araya, R. A. & Harding, A. K. 1999, ApJ, 517, 334 

Araya-Gochez, R. A. & Harding, A. K. 2000, ApJ, 544, 1067 

Avrett, E. H. 1965, SAO Special Report, 174, 101 

Avrett, E. H. & Hummer, D. G. 1965, MNRAS, 130, 295 

Basko, M. M. & Sunyaev, R. A. 1976, MNRAS, 175, 395 

Beloborodov, A. M. 2002, ApJ, 566, L85 

Burderi, L., Di Salvo, T, Robba, N. R., La Barbera, A., & Guainazzi, M. 2000, ApJ, 530, 429 

Caballero, I., et al. 2013, ApJ, 764, L23 

Caballero, I. & Wilms, J. 2012, Mem. Soc. Astron. Italiana, 83, 230 

Chandrasekhar, S. 1960, Radiative transfer (New York: Dover) 

Coburn, W., Heindl, W. A., Rothschild, R. E., et al. 2002, ApJ, 580, 394 

Filippova, E. V., Tsygankov, S. S., Lutovinov, A. A., & Sunyaev, R. A. 2005, Astron. Lett., 31, 729 

Harding, A. K. & Lai, D. 2006, Reports on Progress in Physics, 69, 2631 

Ivanov, V. V. 1973, Transfer of radiation in spectral lines (Washington: US Department of Commerce, National Bureau of Standards) 

Kaminker, A. D., Fedorenko, V. N., & Tsygan, A. I. 1976, Soviet Ast., 20, 436 

Klochkov, D., et al. 2012, A&A, 542, L28 

Kreykenbohm, I., et al. 2008, A&A, 492, 511 

Lyubarskii, Y. E. & Sunyaev, R. A. 1988, Soviet Astronomy Letters, 14, 390 

Mihara, T„ Makishima, K., & Nagase, F. 1998, Adv. Sp. Res., 22, 987 

Misner, C. W., Thorne, K. S., & Wheeler, J. A. 1973, Gravitation (San Francisco: W.H. Freeman and Co.) 

Mitrofanov, I. G. & Tsygan, A. I. 1978, A&A, 70, 133 

Pavlov, G. G., Mitrofanov, I. G., & Shibanov, I. A. 1980, Ap&SS, 73, 63 

Pechenick, K. R., Ftaclas, C, & Cohen, J. M. 1983, ApJ, 274, 846 

Poutanen, J. & Beloborodov, A. M. 2006, MNRAS, 373, 836 

Poutanen, J. & Gierlihski, M. 2003, MNRAS, 343, 1301 

Rybicki, G. B. & Lightman, A. P. 1979, Radiative processes in astrophysics (New York: Wiley-Interscience) 



A reflection model for the cyclotron lines in X-ray pulsars 

Sikora, M., Madejski, G., Moderski, R., & Poutanen, J. 1997, ApJ, 484, 108 

Sobolev, V. V. 1963, A treatise on radiative transfer (Princeton: Van Nostrand) 

Staubert, R., Shakura, N. I., Postnov, K., et al. 2007, A&A, 465, L25 

Suleimanov, V., Poutanen, J., Revnivtsev, M., & Werner, K. 2011, ApJ, 742, 122 

Suleimanov, V. R, Pavlov, G. G., & Werner, K. 2012, ApJ, 751, 15 

Tsygankov, S. S., Lutovinov, A. A., Churazov, E. M., & Sunyaev, R. A. 2006, MNRAS, 371, 19 

Tsygankov, S. S., Lutovinov, A. A., & Serber, A. V. 2010, MNRAS, 401, 1628 

Wasserman, I. & Salpeter, E. 1980, ApJ, 241, 1107 



